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We have developed a microscopic theory of ac Josephson effect in short 
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ballistic superconducting constrictions with arbitrary electron transparency, 
and constrictions with diffusive electron transport. As applications of the 
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g I theory we study smearing of the subgap current singularities by pair-breaking 

qq | effects in the superconducting electrodes of the constriction, and the structure 
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with the proximity-induced gap in the normal layer. 
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It is quite surprising that despite the long history of Josephson junctions, the microscopic 
theory of the ac Josephson effect in junctions with arbitrary transparency is still absent. It 
is known that the basic mechanism of electron transport in such junctions at finite bias 
voltages is the process of multiple Andreev reflections (MAR) jl], but the first attempts 
at quantitative description of MAR based on the Green's function method 0-§J did not 
produce final results. Final progress in this direction has been made recently |||7j but only 
for the case of fully transparent junctions. Quantitative description of MAR in junctions 
with arbitrary transparency is possible within the framework of the Bogolyubov-de Gennes 



(BdG) equations [p|-4tT11 , but this approach does not allow to incorporate inelastic scattering 
and pair-breaking effects in the junction electrodes. This circumstance is crucial because 
dynamics of MAR makes finite- voltage electron transport in high-transparency junctions 
very sensitive to such scattering processes. The aim of the present work was to develop the 
theory of MAR in junctions with arbitrary transparency between the superconductors with 
general microscopic structure. The theory allows us to study, for example, the smearing 
of the MAR-related current singularities by pair-breaking scattering, and the ac Josephson 
effect between the normal conductors with the proximity-induced superconducting order 
parameter. 

The basic model of a high-transparency Josephson junction is a short superconducting 
constriction (shorter than the coherence length £ and elastic and inelastic scattering lengths 
of the electrodes) with a transparency D. The general approach to the description of such 
a constriction requires that we use the Green's function technique. In this technique, the 
constriction is described with the non-equilibrium quasiclassical Green's function G which 
is 4x4 matrix consisting of 2x2 retarded, advanced, and Keldysh matrixes G R ' A and G 
12], ^| . To calculate the current we need to know only the asymmetric part J of the Green's 
function J = G(p Z j) — G(—p Z j). Here p z is momentum in the transport direction, and j = 1, 2 
numbers the constriction electrodes. Solving the quasiclassical equations for G inside the two 
electrodes and matching the solution across the constriction with the boundary conditions 
Q one can show that the matrix J(t, t') is continuous inside the constriction and is given 



by the expression ||: 

J = 2DG_ * G + * (I-DGL * G-)- 1 = G_ * (G + + 2AI)- 1 - (G+ + zAl)" 1 * GL . (1) 

Here the product denoted by * means the convolution with respect to the internal time 
variable, i.e., G^ * G v = J dtiG fl (t,ti)G T1 (t 1 ,t'), and we used the following notations: A = 
y/(l-D)/D, l=l(t-t'), G± = (G 1 ±G 2 )/2, Gjit^') = - t')S*(t'). In the 

last equation, gj(t) = J gj(e) exp(— iet)de/(2ir) is the equilibrium Green's function of jth 
superconductor: &(e) = \gf{e) -gf{e)) tanh(e/2T), gf{e) = gf{e)f z + /f (e)if v = -[gf{e)}*, 
the matrix Sj is : S^(i) = exp[iipj(t)f z /2] where (pj(t) is the phase of the order parameter in 
the jth electrode, and the phase difference if = ip 2 — (fi is determined by the applied voltage 
V(t): <P (t) = (2e/h)V(t). Equation ([1]) gives the following expressions for the Keldysh 
component of J: 

j=(G R *Q 1 -Q 1 *G A -G R *Q 2 + Q 2 * G A )/2 (2) 
= J R *f+-f+*J A + (Gf *P 1 -P 1 *G 2 A -G R *P 2 + P 2 * Gf)/2 (3) 

where Qj = Q R * Gj * Q A , Q» = (G"j + iAi)" 1 , > = G M * & - Q" * Gt, Pj = Q R * 
(Gj */_-/_* Gf) * Q A , f ± = (/i ± / 2 )/2, and f ) = Sj(t)f(t - t')S*(t>) is the matrix 
distribution function of the jth superconductor with /(r) = / tanh(e/2T) exp(— ier)de/ (2tt). 
The current in the constriction can be found from eq. (fj) or (|3|): 

/(*) = ^Trf, J(M) , (4) 

where i?o — vr/i/e 2 for a single-mode constriction. For a constriction with large number of 
propagating electron modes and cross-section area A, angular averaging over momentum 
direction should be taken into account in (|), and R = (2n/epFi) 2 h 3 / A, where ppi = 

Equations (fj), (^) can be used to find the current not only in the ballistic constrictions 
but also in short diffusive constrictions, i.e. channels with large number of propagating 
modes and length d that satisfies the condition I <C d <C £, where I is elastic scattering 
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length. Indeed, in this case solution of the quasiclassical equations for Green's function 
gives for the matrix J the expression J = (2lvp z j dvp) hi(G 2 * G\) 0, which reduces to the 
following form [13[ : 

Alv Fz f°° d\ 



J 



After the substitution: A = J (1 — D)/D, we get for the current 



where J(t,t;D) is given by eq. (|2|), and Rn is the normal-state resistance of the channel. 



Equation (jg) shows that similarly to the approach based on the BdG equations JTTJ, in the 
general Green's function approach, the current in the diffusive superconducting constriction 
can be written as a sum of independent contributions from infinite number of ballistic 



propagating modes with the distribution of transparencies given by the Dorokhov's WA 



density function (nTi/2e 2 R N ) / D^Jl — D. 

To calculate the current at arbitrary bias voltages across the constriction we should 
transform eqs. @) and (§) further. Separating the normal-state contribution V/Rn, where 
Rn = Ro/D is the normal-state junction resistance, we obtain the following symmetrized 
expression for the current: 

I(t) = ^- + 6I 12 (t) + 6I 21 (t), (7) 

SI jk (t) = -^b[UG? * Qjk - q jk * Gf))(t, t) , (8) 

where q jk = qf k * g k * qf k , f jk = 2/ (G^ + <% + 2iXi)~ 1 , and g k = gj? * f - f * g£. In equation 
(|), all distribution functions correspond to zero potential, and in contrast to the previous 
definition, the function G^ and G 2 are defined similarly, Gj = S(t)gj(t — t')S*(t'), where 
S(t) = exp[i<p(t)f z /2]. Different representations for the current can be obtained from eq. 
(|). We will use the representation in which the Fourier components qf k ,A (e, e') are expressed, 
respectively, through the matrices ocf k (e, e') and af k (e,e') = [af k (e', e)]t, where, for example, 
af k obeys the following equation: 



(Pf-f f * r/j?) * &f k = ri+ff , (9) 



In this equation, r = yl — D is the reflection amplitude of the constriction, and other 
notations are: p%{t,t) = i - rTf (t - t')r x , i)j?(t,t') = -rl + if (t - t')r x , T%(t) = 
J^(e)eM-tet)de/(2n), f£(M') = T*(t - t') exp{zf z [^) + <p(t')}/2}r x , >y*( £ ) = 
fk'( e )/[9k( e ) + 1]> 7fcK e ) = 1k*( e )- As we wm see later, j k (e) has the meaning of amplitude 
of Andreev reflection at the interface between the constriction and the kth superconductor. 
In terms of ajjj. eq. (|8|) reduces to the following form: 

MM = 7^r f^:f deJ jk (e, e - U )e~^ , (10) 



Sei^n J 2tt 

where 



J jk (e,e') = Trir y [^(e)W k (e')af k (e,e f ) -^(^O] 
dei 



7 fc R (e) 2 



Equations (|T0|) and (0) show that the problem of finding the current in short ballistic super- 
conducting constriction for arbitrary time-dependent bias voltage reduces to the problem 
of solving eq. (^]), which is a Fredgholm integral equation (see, e.g., [|15|). Similar pro- 
cedure applied to (|6|) gives the current in diffusive constriction. It should be noted that 
this approach assumes that all frequencies (frequency of Josephson oscillations and typical 
frequency of voltage variations) are much smaller than the inverse of the time of electron 
motion through the constriction. 

Equation (H) can be solved easily in the case of the dc bias voltage V, when the phase 
difference is: (p(t) = ujjt + <p , where ujj = 2eV/h is the Josephson oscillation frequency. In 
this case solution of eq. (|§D can be written as a series , <x,-fc(e, e') = 27iJ2n^n(jk)(^)^(^ — 
e' — nhujj) (we omit the superscript R), in which the amplitudes a n (jk) are determined by a 
system of recurrence relations. To obtain these relations explicitly it is convenient to write 
the matrix &j k as 
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( b + o- ^ 



i a + b j 



m 

Equation @ shows that the pairs of functions dj k , bj k with s — ± satisfy the same equations 
but corresponding to different polarity of the bias voltage. For dc voltage we get the following 
recurrence relations for the matrix elements of a n fjky. 

a n +i-lj^2n+i)lk^2n)an + r[y j{e 2n+ i)b n - ^k^2n+2)b n +l\ = 7i( e i)<W (11) 
c(e 2n +i)b n +i - d(e 2n )b n + c(e 2n -i)6„-i = -r5 n0 , 

where we used the simplified notations a n = a^ k ^(e) and b n = b^(j k Je), and other notations 
are: e n = e + neV, c(e) = D lj {e) lk {e + eV)/[l - j* (e)], d(e) = 1-^ (e) +D 1 ](e + eV) /[l - 
7?(e + eV)] + Z>y|(e)[l — 7? (e — eV)]. The amplitudes a„y/») and & n (jfc) determine Fourier 
components of the current I(t) = J2n^nGxp(inujt), which according to eqs. (|7|) and (|TUp are 
given by the expression J n = VS n>0 /R N + I n (i2) + /n(2i), where 

1 f 00 

W) = 15 77" / deW k (e)$ n{jk) (e), (12) 
zeitQ J -00 



El 



'n(jfc)(e) = -7fc(e2n)an(jfc)(e) - 7fc(e„ 2n )al n(jfc) (e) + 



1 + 7fc(e2m+2n)7fc( e 2 m )][&n+m&m ~ a n+m a^J (jfe) (e) , 



For a constriction between two BCS superconductors the recurrence relations (|TT|) and equa- 
tion fll2|) for the current reduce to the corresponding expressions that can be obtained from 
the BdG equations - see ||, where these expressions were derived for symmetric constric- 
tion. This means that for the purpose of description of the ac Josephson effect in a short 
constriction all information about the microscopic structure of a superconducting electrode 
of the constriction is contained in the function 7 R (e) (introduced after eq. (|9|)) which has 
the meaning of the amplitude of Andreev reflection from this superconductor. 

Equations QTTJ) can be solved by standard methods |T5| (see also [§]). Namely, it fol- 
lows from flTlD that b n (e) = bo(e) Ilm=±i P±( e 2m), where the functions p±(e) are solutions of 
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the equations p±(e) = c(e T i)/[d(e) — c(e±i)p±(e±2)]; the two signs (±) here correspond to 
n > 1 (n < —1), and 6 ( e ) = r [^( e ) — c ( e i)P+( e ) ~~ c ( e -i)P-( e )] _1 - These relations provide 
the basis for convenient numerical evaluation of the current. They also can be used to find 
analytical solutions in certain cases. In particular, in the low voltage limit, eV <C Aj, we 



get explicitly p±(e) = re(e) - *//e a (e) - 1 where «(e) = {[l- 7l 2 (e)] [1 - 7 2 (e)] + Djf(e) + 



D 7 2 ( e )}/2D 7l (e) 72 (e), and fc (e) = r/2L> 7l (e) 72 (e)^K 2 (e) - 1. Thus, eqs. © and (pD 
enable us to find the current in the constriction between the superconductors with arbi- 
trary quasiparticle spectrum. As was discussed above, this feature of our approach is very 
important, since all deviations of the quasiparticle spectrum from its "ideal" BCS form af- 
fect strongly the subharmonic gap structure and especially the low-voltage behavior of the 
current. 

The quasiparticle spectrum of a superconductor can differ significantly from the BCS 
form, in particular, due to pair-breaking effects. One of the most important example of these 
effects is scattering on paramagnetic impurities - see, e.g., P^JT7| . As a first application 



of our general theory we consider the ac Josephson effect in a constriction between two 
superconductors with paramagnetic impurities. The retarded Green's functions of such 
superconductors are determined by the relations pi : 



R 
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- = u[l-^==], (13) 



u 



where ( is the pair-breaking parameter, ( = h/r s A, with t s being the spin-flip scattering 
time. Similar expressions for the Green's functions are valid also in the situations with 
other pair-breaking effects for instance, for dirty thin superconducting film of thickness 



ds £ in magnetic field H parallel to the film. The Green's functions of such a film 
are given by eq. ( |1"3"D with the pair-breaking parameter ( = IvpieHds) 2 /18A. Weak pair- 
breaking effects result in smearing of the BCS singularity in the quasiparticle spectrum and 
suppression of the superconducting energy gap to a reduced value A g = A(l — £ 2 / 3 ) 3 / 2 . The 
gap disappears completely at £ > 1. Figure 1 shows how these changes in the quasiparticle 
spectrum are reflected in the current- volt age (I — V) characteristic of the constriction. All 



gap-related features in the I — V characteristic are rapidly broadened at small (, and it 
becomes practically linear in the gapless regime £ > 1. 

As another application of our theory we consider a constriction between two normal 
conductors in which superconductivity is induced by the proximity effect, i.e. an S/NcN/S 
junction. Besides general interest to the proximity effect, the importance of this model is due 



to its relevance for realistic description of the high critical current tunnel junctions [ lq| . We 
consider a particular case of a thin dirty iV layer of thickness d n -C £ n with low transparency 
of the S/N interface (D 1 ) < 1. It is assumed that the resistance of the interface is still 
negligible in comparison to the constriction resistance. [|19| The Green's functions in the N 
layer are given then by the first equation in (|T3|) with u = (t+ilbQs) l^lbfs ~~ see > e -S-' IPIPll 
and references therein. Here g§ and f$ are the Green's functions of the superconductor, 
and 7fe//i = (D') vp n /kd n is the characteristic tunneling rate across the S/N interface which 
is assumed to be larger than electron-phonon inelastic scattering rate in the iV layer. The 
energy gap A 9 is induced in the N layer due to the proximity effect. If the 5* electrode of the 
structure is the BCS superconductor with energy gap A, the induced gap A 5 is determined 



by the equation A g = A7 b /[y A 2 — A 2 + 7 b ]. Existence of the induced gap implies that there 
are two peaks in the density of states of the iV layer, at energies A g and A. This structure 
of the density of states results in a complex structure of the subharmonic gap singularities in 
the I — V characteristics of the constriction. Example of such a structure is shown in Fig. 2 
for 7b/ A = 1, when A g ~ 0.54 A. We see that the most pronounced current singularities in 
this case are the subharmonic singularities at eV = 2A g /n associated with the induced gap. 
Also visible are the singularities at "combination" energies A + A g and (A + A g )/2. 

In the case of fully transparent constriction with D = 1 we can obtain explicit analytical 
expressions for the current since the recurrence relations ( [TT]) can be solved directly in 
this case. In particular, in the low voltage limit eV A g f(A g ) we get for symmetric 
constrictions: 



w 2eR N J u 2 {e) - cos 2 (v?/2) ' v ; 



F ± {e,V)=[F(e,V)±F(e,-V)]/2, 



F(e , ± v) = M - £ a/ 4£> ex P (- J; ' § m \^F) 

These expressions imply that the current can be written as a sum of two contributions, one 
from the two discrete states inside the energy gap with energies e± = ie^, determined by 
the equation ^(e^) = |cos(<^/2)|, and another contribution from the continuum of states 
above the gap. In contrast to constrictions between the BCS superconductors |22].[23].|7| , in 
the general case of the non-BCS spectrum of quasiparticles both the discrete and continuum 
contributions to the current are significant. In particular, for the ballistic constriction in 
presence of pair-breaking effects one finds from eqs. (|T3|) and ([14]) that the current Id carried 
by the subgap states is: I d = I |sin(<^/2)| (l — (/ sin 3 (<^/2) J 6 (jsin(<^/2)| — C 1 ^ 3 ), where 
/ = 2sgn(V)Agf(A g )/eRN- At ( = this expression reduces to the one found in 0. The 
total current I in the constriction as a function of the (time-dependent) phase difference p 
calculated from eqs. (|13|) and (|TJ]) is shown in Fig. 3. We see that even relatively small ( 
has strong effect on the dynamic current-phase relation, suppressing the dc component of 
the current and making it more similar to the stationary current-phase relation. 

In conclusion, we have developed the microscopic approach to the calculation of current 
in short ballistic and diffusive constrictions between the superconductors with arbitrary 
quasiparticle spectrum. We used the developed approach to study the time- dependent cur- 
rent and the dc current in constrictions between superconductors with the pair-breaking 
effects, and also between normal conductors with the proximity- induced superconductivity. 

This work was supported in part by the ONR grant # N00014-95- 1-0762, Russian Fund 
for Basic Research, grant No. 96-02-18613, and by the U.S. Civilian Research and Develop- 
ment Foundation under Award No. RP 1-165. 
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FIGURES 



DC current-voltage characteristics of a short ballistic constriction with transparency 
D between the two superconductors with pair-breaking effects. The curves are shifted for 
clarity along the current axis and illustrate the smearing of the subharmonic gap structure 
with increasing strength of the pair-breaking. The values of the pair-breaking parameter ( 
are (from bottom to top): ( = 0.01, 0.1, 0.3, 1.0. The upper curve with ( = 1.0 corresponds 
to the regime of gapless superconductivity. 

DC current- voltage characteristics of a short symmetric S/NcN/S constriction for dif- 
ferent values of the constriction transparency D. The curves show the complex subharmonic 
gap structure associated with the two energy gaps: A in the S region, and proximity-induced 
gap A g in the N region. Parameter 7^ characterizes the transparency of the S/N interfaces. 
For discussion see text. 

Dynamic current-phase relation of a short symmetric constriction between supercon- 
ductors with pair-breaking effects at low bias voltages, vanishing temperature, and different 
values of the pair-breaking parameter. 
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